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The three-dimensional Gross-Neveu model in R 1 x S 2 spacetime is considered 

at finite particles number density. We evaluate an effective potential of the 
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R and nonzero chemical potential /i. We then derive the critical values of 

(R, jj) at which the system undergoes the first order phase transition from 

the phase with broken chiral invariance to the symmetric phase. 
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I. INTRODUCTION 



Last time four-fermion field theories in (2 + l)-dimensional Minkowski spacetime, which 
are known as Gross - Neveu (GN) models [[I], are under extensive investigation for purely 
theoretical motivation and also due to their applications to planar condensed matter physics. 
Such theories possess many desirable properties: the renormalizability in 1/N expansion, 
dynamical breaking of chiral symmetry and generation of fermion mass for a large coupling 
constant as in QCD M, the analogy to the BCS theory of superconductivity in two spa- 
tial dimensions and the possibility to describe the new phenomenon of high temperature 
superconductivity |§, the reduction to the S = 1/2 quantum antiferromagnet Heisenberg 
model in the continuum limit [§J and so on. Main features of these models, obtained in large 
N expansion technique, are confirmed within the framework of other nonperturbative ap- 
proarches |J. Since there are no closed physical systems in nature, the influence of different 
external factors on the vacuum of the simplest GN model was considered. In || some critical 
phenomena of this theory were studied at nonzero temperature T and chemical potential 
/i. Recently, on the same foundation a new property of external (cromo-)magnetic fild H 
to promote the dynamical chiral symmetry breaking was discovered [0]. (At present it is 
the well known effect of dynamical chiral symmetry breaking catalyst by external magnetic 
field ||, which is under intensive consideration ||.) The role of T, H as well as of //, H 
in the formation of a ground state of the GN model was also clarified ||10|| . The study of 
dynamical symmetry breaking in spacetimes with curvature and nontrivial topology is also 
of great importance, since in the early universe the gravity was sufficiently strong and one 
should take it into account. There is a reach literature on this subject (see the review [pi ). 



The papers in [|T2| , p!3| are the first ones where the effect of curvature and nontrivial topology 
on the chiral symmetry breaking in four - fermion models was discussed. The curvature - 
induced first order phase transition from a chiral symmetric to a chiral nonsymmetric phase 
was shown to exist in that models in the linear curvature approximation. It turns out that 



in specific spacetimes such as Einstein universe jnj and maximally symmetric spacetimes 
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15] the above mentioned models may be solved exactly in the leading order of large N 
expansion technique. Finally, dynamical symmetry breaking in the external gravitational 
and magnetic fields was considered fTEf . It is well - known that low dimensional four - 
fermion field theories, especially the (2+1) - dimensional GN model, in curved spacetimes 
P~3| , [T^JTT| , P~5| ] and in the nonsimply connected spacetimes [H||2T)| | may be very useful for the 



investigation of physical processes in thin films and in the materials with layer structure. 
The matter is that external stress, applied to the planar system, may change topology and 
curvature of a surface. The consideration of above mentioned low dimensional models in 
curved spaces is an interesting subject by itself, since in these cases the quantum theories 
can be computed exactly. The great amount of observable physical phenomena are due to 
nonzero particle density (superconductivity, quantum Hall effect etc.). So in the present 
paper the influence of both chemical potential and curvature of space on the phase structure 
of (2+1) - dimensional GN model is studied. Especially, we shall consider R 1 x S 2 space- 
time to clarify our discussion (the effect of temperature on the vacuum of GN model in this 
spacetime was considered in ||TT|). In Section II, we evaluate the one-loop effective potential 
in R l x S 2 spacetime at nonzero chemical potential. Presented in Section III is the detailed 
analysis of the effective potential, which shows the existance of a phase transition restoring 
the chiral symmetry of the system as the curvature R and chemical potential /i are varied. 
Finally, we summarize our results in Section IV. 

II. EFFECTIVE POTENTIAL IN R 1 x S 2 SPACETIME AT (i / 

The R 1 x S 2 space-time is chosen to be 

ds 2 = dt 2 - a 2 (d6 2 + sin 2 9 d<j> 2 ), (1) 

where a is its constant radius. The four-fermion model in this spacetime is described by the 
action |TT|JT^| 

A 2 



S = / d 6 x^g 



#^(*)v^ + ^(^) 2 



(2) 



where V M is the covariant derivative and the summation over j is implied (j = 1,2, ..,N). 
Here fermion fields ipj are taken in the reducible four dimensional representation of SL(2, C). 
For this case the algebra of the 7 - matrices is presented in . This action has the discrete 
chiral symmetry, 

tj; -> 75^. (3) 

As a result, the chiral symmetry is maintained at any order of ordinary perturbation theory. 
However as can be seen in different nonperturbative approaches |],|2||| the symmetry may be 
broken dynamically for large values of coupling constant A. To see the nonperturbative fea- 
tures such as spontaneous symmetry breaking and dynamical mass generation in the present 
model, it is convenient to rewrite above action in an equivalent form |J by introducing the 
auxiliary field a(x), 

N 



S = Id Xy/—g 



(4) 



2A 

This expression explicitly says that the vacuum expectation value of a field plays the role 
of mass for the fermions. In order to find the effective potential in the theory with the 
action Eq. (2) we follow p2]fT3| where this quantity was considered in a weak curvature 



approximation. First of all let us integrate over the fermion fields in Eq. (4) and evaluate 
an effective action S e ff(a) describing the self-interaction of a field: 

exp(iNS eff (a)) = J DipDipexp[iS(ip,ij;,(T)]. (5) 

Here we use the 1/N expansion which is the fermion- loop expansion. In the mean-field 
approximation, where the cr(x) field is assumed to be constant, and to the leading order in 
the large N, one can obtain the one-loop effective potential U(a) from the action S e ff(a): 

2 

U(a) = ^+i tr{x\ ln(z 7 ^(a;)V M - a) \x), (6) 

where tr is over indices other than spacetime indices. Using the Green function Gp(x,y; a) 
defined by the relation 



G F {x,y,(r) = {x\{i^V lt -(T)- x \y), (7) 
we rewrite Eq. (6) as follows: 



a 2 



U(a) = — —i tv\nG F (x,x;a). (8) 
The logarithm may be eliminated from this equation by introducing the parameter s: 



In 



K-a- 




[ K 


= -/ 

•/ 



ds-^L., (9) 



where an operator K is given as i7 M (x)V Al in the present case. Therefore, Eq. (8) is rewritten 
in the following form: 

where the momentum-space Green function G F (k; s) has been used. Now one may introduce 



the Riemann normal coordinate [21] with origin at any point in the spacetime. In this 



coordinate system we use the following expression for the Green function G F (k; s): 

Lrp{K,S)— -,rw,„2 „2N2 + o U P» K K 



k 2 -s 2 12 {k 2 — s 2 ) 2 3 ^ (k 2 -s 2 f 

~lp a J Cd Rcda^ _ ^ 2 , (11) 

where J ab = |[7 a ,7 fc ], and the Latin and Greek indices refer to a local orthonormal frame 
and general coordinate system, respectively. Eq. (11) is the linear approximation for the 
Green function G F (k;s) in the curvature R |TT|--[l3l. According to the well known method 
developed in p2| |, to obtain eq. (11) one should neglect any terms involving derivatives 



higher than that of the second order in the metric tensor expansion. Now let us consider 
the effect of nonzero chemical potential \x on the system. As is well - known, the fermion- 
number density is directly related to the chemical potential /i. Mathematically, the presence 
of nonzero chemical potential is realised by shifting the energy levels ko in the propagator 
G F (k;s) by the amount of fi Thus, we are in order to evaluate the effective potential 



U(a) in Eq. (10) under effects of both R and /i. Using the contour integration method 



23| , we can perform the integration over momentum k^. Denote ii the integral of the first 
term in Gp(k;s) over k and s. Its calculation proceeds as follows: first, the procedure of 
integration over k , denoted as /{, gives the result: 

T'd \ + f dk j°(k + fM) +fh + s 
I^k, s) = tr ' 



(12) 



2vr (k + /i) 2 - El 

2 r io ° s dz 2 r s dz 
+ - 



TT J— ioo Z 2 — Ell 



El 



TT JC Z" 

2 r°° s dz 2is „, 
E k E k 



TT J —ioo 

Here, El = k 2 + s 2 , the contour C is given in Fig. 1, and the unit step function 8(x) 
1 for x > 0, 9(x) = for x < has been used. Thus we get 

d 2 k 



I\ = —% I ds 
Jo 



(2tt 



a 



a_ A_ 

3/T 7T 2 



+ - a) 



V 2 

2^ 



1 

3tt' 



-<7 



(13) 



where A is the cutoff parameter. Here and in the following discussions, we may confine 
ourselves to the o > region due to a reflection symmetry a «-> — er of the effective potential 
U(a). However, note that this symmetry is broken when the system selects one of the two 
ground states. In a similar way one finds the contributions 7 2 ? 1 a of the remaining terms 
of Gp(k; s) to the potential (10): 



h 
h 



R 

24^ 
1 

127T 

1 

~ 12V 
h = 0. 



-a + 6(ji - a) (a - ^a 2 ) + 9{a - 



(14) 



In the third line of Eq.(14), we have used the relation R 0Q = 0, which follows from the metric 
Eq.(l) of the spacetime under consideration. However, the fourth line of Eq.(14) is due to 
a relation trpy'T^ ] = 0. At this stage it is convenient to introduce the mass parameter M 
instead of the coupling constant A by the following way H]: 
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1 

A 2 



r A d 3 k E 1 
J (2vr) 3 k\ + M 2 




(15) 



So, we shall consider the case A > A c only, where A~ 2 = 4 J A (i 3 fcB(27r)~ 3 fc^ 2 . Summing up all 
terms in Eq. (13) and (14) and inserting above equation into Eq. (10), one sees that the 
two A-dependent terms cancel out, and thus the finite effective potential to one- loop order 
is obtained. Then, the \i- and independent one- loop contributions Uj^Ja) to the potential 
U(a) are completely separated from the Minkowski-space result: 



where Up(cr) is the effective potential of the original theory in flat Minkowski spacetime. 
Here 



In this expression one may find the following two facts. Firstly, Uj> u (a) is finite and, as 
R, [L — > 0, UftJa) vanishes. Thus the renormalisation procedure is identical to the case of 
Minkowski spacetime. Secondly, in the limit R — *■ 0, U(a) is reduced to the Minkowski- 
space effective potential Uf{&)- It is well established that there are two distinct phases in the 
three - dimensional GN model (l|,@,0. For a weak coupling phase with the coupling A < A c 
we have (a) = 0. Thus the fermions are massless and the chiral symmetry remains intact. 
However, for the strong coupling phase A > A c , a field has nonzero vacuum expectation value 
(a) = M , so the chiral symmetry Eq. (2) is dynamically broken and fermions aquire mass, 
which is equal to the mass parameter M from Eq. (15). For the simplicity of our analysis in 
the next sections, we shall introduce the following rescaled dimensionless quantities defined 
as U(x) = 7iU(cr)/fi 3 , R = R/fi 2 , x = er//i, and p, = \ijM. In terms of these quantities, Eq. 
(16) is rewritten in the much simpler form: 
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U(a)=U F (a) + UUa) 



(16) 




+ —e(a-fi) M/i 2 + -) 



(17) 



'1 _ I _ + Bs. for r < 1 

£/(*) = < A ^2_+i2> (ig) 



3 )T + f + K 1 + f) ) forx>l, 



where one sees that U(x) is a continuous function at a; = 1. We also wish to find the 
induced fermion mass (a) as a function of curvature R and chemical potential \x. Then, the 
gap equation for the fermion mass can be obtained by taking the derivative of the effective 
potential U (x) with respect to x, and so we obtain 







(! - 7 ~ £) x + for x < 1 

(19) 

3:2 ~~ f + H' for x > 1. 



III. RESTORATION OF CHIRAL SYMMETRY 

Now we shall analyze in detail the effective potential Eq. (18) in order to investigate 
the phase structure of the model in the (R, /z) plane. The fermion mass (<r) will be derived 
which depends on R and fi and the nature of the phase transitions will be discussed. To 
clarify our discussion, we shall consider three distinct cases: ji ^ and R = 0, then R ^ 
and fi = 0, and finally R ^ and /i ^ 0. 

A. The case fi / and i? = 0. 

Let us first examine the effect of nonzero chemical potential on the system. In the limit 
R — > 0, the effective potential Eq. (18) is reduced to a simple form: 

1-4)^, for x < 1 



C/(x) = <{ v 2 (20) 

§)^ + |, forx>l. 

To see a phase transition as Jx increases from a broken phase to a symmetric one, it is 
necessary to examine the behavior of U(x) as a function of Ji. It is possible to find the 
following two properties of U(x). For p, > 1, U(x) is a monotonically increasing function of 
x, and so the global minimum of U (x) occurs at x — 0. While for // < 1 U (x) has a global 
minimum at x — 1/jx with the value: 



These facts indicate that the system undergoes a phase transition from the (a) = M state 
to the (a) = state at the critical value fi c of the chemical potential, given as 

H C = M. (22) 

Solving the gap equation for the induced fermion mass, Eq. (19) with R = 0, one can find 
that 

(a) = M (23) 

below /i c , and (a) = above \i c . Except at ii — /i c , the order parameter (a) does not depend 
on the value of fi. That is, the value of order parameter (a), which minimizes the potential, 
jumps discontinuously from a = M to a = at the transition point /i c . Hence, at the point 
H — fi c we have a first order phase transition from a massive chirally broken phase to a 
massless chirally invariant phase of the model. 

B. The case R ^ and = 

In this case only the effect of curvature on the system will be considered. In the limit 
/i — > 0, the general effective potential Eq. (18) has the following form: 

U(a) = U F (a) + U\{a) 



a 



2 / 3 \ R 



— \a--M + a. (24) 

3tt V 2 / 24tt v ; 



This expression coincides with that obtained in |TJ| . ^From Eq. (24) one can see that in the 
region of small values of a the dominant contribution to U(a) comes from the /^-dependent 
linear term in a. Thus, there is a potential barrier between a = and second local minimum 
of U(a). As a result, it turns out that as the curvature R increases the discontinuous phase 
transition occurs from a chirally broken phase to a symmetric one. The critical value of the 
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curvature R c , at which a first order phase transition occures, is determined by the following 
two conditions: 



where a denotes second nonzero local minimum of the potential. Furthermore, one may 
find that only for R > R c the minimum of the potential at the symmetric point a = is 
lower than asymmetric local minimum at a nonzero &q. ^From the gap equation Eq. (19) 
with fj, — one can evaluate the local minimum of the potential ctq, 



which at the same time equals to the fermion mass (a), induced under the influence of 
curvature R for R < R c only. Thus, applying the critical condition Eq. (25) to the effective 
potential Eq. (24), one can obtain the critical curvature 



The phase transition under the influence of R is a first-order one since it occurs discontinu- 
ously. 

C. The case R ^ and fj, ^ 

In given Subsection we are going to explore the general case when the system is specified 
by the curvature and finite chemical potential. To investigate the vacuum structure of the 
system as R and \i are varied, one must first examine the behavior of the potential U(x) 
as a function of R and /i. It is very helpful to sketch qualitatively the effective potential 
U(x) from Eq. (18). For p > 1 (fi > M) the global minimum of U(x) occurs only at x — 0. 
While for p < 1 (/i < M) , the global minimum of U (x) lies at nonzero point certainly. 
Therefore, when p, < 1, it turns out that the system undergoes a phase transition from the 
(a) 7^ vacuum state to the (a) = state at certain critical curvature R c depending on 
\i . Using a much detailed analysis of the effective potential U(x) in Eq. (18), one can 



U'(a ) = and U(a ) = 0, 



(25) 




(26) 



R c = 4.5 M 2 . 



(27) 
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see that until the system approaches the critical point with the increase of curvature, the 
second local minimum of U(x) occurs only in the region x > 1. Therefore, in the procedure 
of determining the critical value of the curvature R c , the effective potential needs to be 
considered only in the x > 1 region in Eq. (18). 

In this case we can obtain the critical curvature R c also using the condition given in 
Eq. (25), with the only change <jq — > x , where Xq denotes the second local minimum of 
the potential. That is, in the present case the phase transition under investigation is also 
a first-order one. As can be easily checked from the gap equation Eq. (19), the second 
minimum lies at the point 



^H 1- ^'- (28) 

Thus, the critical condition Eq. (25) with this value for x leads to the self-consistent relation 
on the critical curvature R c : 

24r 2 

16x1 — + ( 2x o + l)Rc + 8 = 0, (29) 

where x has the value given in Eq. (28), with R replaced by R c . The numerical solutions 
of Eq. (29) are illustrated in Fig. 2. Note that as /i — > the R c approaches 4.5 M 2 and 
as R — > 0, the ji c approaches M. These limiting cases have been already discussed in the 
previous Subsections. Eq. (28) says that the induced fermion mass (a), with (a) = /ix , 
does depend on the curvature R only. That is, (a) does not depend on /i, and thus it has 
the same expression as Eq. (26). In Fig. 3, the effective potentials are given for four distinct 
values of R at fixed \i = 4^-. 

IV. SUMMARY 

In the present paper we have derived the effective potential of the three-dimensional 
Gross-Neveu model in the curved spacetime of the form R 1 x S 2 and with taking into 
account the chemical potential /x as well. Then, the critical curvature R c at which dynamical 
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symmetry breaking disappears has been determined in terms of the induced fermion mass 
M in the limit R,fi—>0 and at nonzero chemical potential /i, as given in Fig. 2. As the 
chemical potential \i increases, the fermion-number density increases also. In Subsections 
A and C, it has been shown that the high density influences the symmetry behavior of the 
system, and so at the critical value fi c or at the corresponding critical number density the 
chiral phase transition is occurs. Then, we have observed that the order parameter (cr) of 
the phase transition, corresponding to the minimum of the potential, does not depend on 
the value of /z, except at the critical value /i = /i c , even though the phase transition is 
induced by the chemical potential. This phenomenon is connected with the fact that the 
composed field a ~ ipip is a real field and carries no charge. It was observed also in two - 
dimensional GN model in R l x S 1 spacetime ||20|| , however, in that model there is another 



massive phase, in which fermions mass is fi - dependent quantity. In Subsection B, we have 
shown that R c = 4.5M 2 . This can be roughly seen from the following two facts. Firstly, on 
dimensional grounds the critical curvature R c must be proportional to the square of some 
quantity with the dimension of mass. Secondly, the effective potential for the composite a 
field in Eq. (24) has two parameters R and M, and so the remaining parameter apart from 
R in this theory is M. Note, that our value for R c is valid only in a weak curvature limit, 
and thus its more accurate value can be obtained by considering higher order corrections to 
scalar curvature R. However, in such improved schemes, it is expected that the system still 
show the same qualitative properties as those found in the previous Sections, including the 
occurrence of a first order phase transition. Finally, one may consider the case of negative 
curvature since the present method has the advantage of being applicable to any metric. 
Then, Eq. (26) indicates that under the effect of negative curvature R the minimum of the 
potential moves farther from the origin than without the curvature effect. Therefore, in this 
case the symmetry restoring phase transition does not happen. 

We hope that the above results may be useful for condensed matter physics as well, and 
for astrophisical applications, especially for the desciption of different phenomena in the core 
of neutron stars. 
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Figure Captions 



Fig.l. The contour C in the complex k° plane. Fig. 2. The critical curvature R c /M 2 as 
a function of nonzero chemical potential fi/M. In region B, chiral symmetry is broken and 
fermions acquire dynamical masses, while in S, the symmetry is restored by the curvature 
effect, and fermions become massless. Fig. 3. The effective potential 7rU(a) /M 3 as a function 
of a/M at the fixed value of fx/M = 1/2. Four interesting cases of R, where R = R/M 2 , are 
considered, and the critical curvature R c is then numerically obtained: R c = 2.96. 
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